Schrodinger equation with a spatially and temporally random potential: Effects of 
cross-phase modulation in optical communication. 
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We model the effects of cross-phase modulation in frequency (or wavelength) division multiplexed 
optical communications systems, using a Schrodinger equation with a spatially and temporally ran- 
dom potential. Green's functions for the propagation of light in this system are calculated using 
Feynman path-integral and diagrammatic techniques. This propagation leads to a non-Gaussian 
joint distribution of the input and output optical fields. We use these results to determine the 
amplitude and timing jitter of a signal pulse and to estimate the system capacity in analog commu- 
nication. 
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I. INTRODUCTION 



A. General 



The ability to transmit information is ultimately lim- 
ited by signal distortion. Information theoryM quantifies 
the extent to which such distortions inhibit communica- 
tion; it was originally developed for the study of radio 
communication or electrical communication along cop- 
per wires. In these cases the signal propagation is linear- 
the received signal is a linear function of the transmitted 
signal. Distortion arises due to the addition of extrane- 
ous signal fluctuations, which propagate linearly along- 
side the original signal. These fluctuations may come 
from noisy amplifiers or other circuit elements, or from 
cross-talk with other messages. Such linear systems with 
additive noise are very well characterized from an infor- 
mation theoretical perspective. 

In modern communication systems, the situation is 
rather more complicated. The very high transmission 
rates- particularly in optical communication- require that 
the transmission medium be operated in regimgs where 
the signal propagation is substantially non-lineaitl. These 
nonlinearities can ]£ad 1[ tQ ir a variety of new mechanisms 
of signal distortioroj iMlJt] and are often very difficult 
to characterize analytically. This has limited the under- 
standing of non-linear channels, particularly from an in- 
formation theoretical perspective. 

Faced with these difficulties, there are two ways in 
which progress may be made. One approach is to perform 
detailed numerical simulations of the underlying partial 
differential equation. This has been the approach of much 
of the work in the literature. However, it is very dif- 
ficult to use the results of such analyses to understand 
the information theoretical aspects of communication in 
non-linear media. An alternative approach is to harness 
physical understanding of the signal propagation to moti- 
vate simple phenomenological models that capture some 
aspects of the non-linear propagation, but which are nev- 
ertheless analytically tractable. This approach has been 



used to understand a numher-aflaspects of non-linear 
propagation in optical fibersQEnffl'QEI. Only recently have 
such approaches been applied to ^udy_thc information 
theoretical aspects of the problemE!JiIdli3. 

When one is interested in optimizing the design of a 
communications system, it is often useful to characterize 
the various mechanisms of signal distortion according to 
their physical origin in the system. This is not always 
the best way of classifying signal distortion from an in- 
formation theoretical point of view. A more convenient 
classification is to divide signal distortions into additive 
and multiplicative noise. The signal distortions in opti- 
cal fibers may be classified in this manneniii It is the 
multiplicative noise that presents the main difficulty in 
understanding non-linear channels. It lies behind sonae 
of the classic unsolved problems in information theorya. 
It is to this aspect of the problem that path-integral and 
diagrammatic techniques are particularly suited. In the 
system studied in this paper, multiplicative noise appears 
as a random potential in a Schrodinger equation. 

The tension between additive and multiplicative noise 
is similar to the tension between disorder and interac- 
tions in many-body physics. In the latter case, the con- 
sequence of this tension is that a full analytical charac- 
terization is not possible. Instead, one must resort to 
approximation schemes that are appropriate under dif- 
ferent regimes of system parameters. It is likely that for 
the same reasons one will be forced to use similar ap- 
proximate schemes to study the propagation of signals 
in non-linear media. Physicists have worked hard to de- 
velop such schemes and some of them are directly appli- 
cable to the study of noHrlinear signal propagation; for 
example, Falkovich et ata have used the theory of opti- 
mal fluctuation to understand the fluctuations in soliton 
amplitude and timing that occur in optical communica- 
tion. In this paper, we use,Feynman path-integralo and 
diagrammatic techniqucsr-3 in order to understand some 
of the effects of non-linearity upon optical communica- 
tion. Most of this paper will concentrate upon charac- 
terizing the propagation of light in the fiber. Towards the 
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end of the paper, we will harness this characterization to 
understand some information theoretical aspects of the 
problem. 



B. The wavelength division multiplexed (WDM) 
optical fiber. 

In this paper, we will consider some of the effects of 
non-linearity upon communication using frequency divi- 
sion multiplexed optical fibersa. Optical fibers have an 
enormous transparent bandwidth. It is far beyond cur- 
rent electronics to modulate at these frequencies. In or- 
der to utilize as much of the bandwidth as possible, it 
is divided into a number of non-overlapping frequency 
bands, or sub-bands, each of which is modulated indepen- 
dently. This is called frequency division multiplexing or 
alternatively wavelength division multiplexing(WDM). 

Light propagating in a WDM optical fiber experiences 
several sources of signal distortion and non-linearity. The 
first of these is due to dispersion; the different frequency 
components of a signal travel at different velocities along 
the fiber, leading to a spreading of a signal pulse. This 
spreading is a linear operation, however, and does not 
reduce the capacity to communicate information. The 
effects of dispersion could be taken into account by a suit- 
able linear transformation of the received signal during 
processing. Usually, however, dispersion is accounted for 
by introducing dispersion compensating elements into the 
transmission link; Spans of standard, dispersive fiber are 
interleaved with spans of dispersion compensating fiber 
in which the sign of dispersion has been reversed, so that 
more slowly moving frequency components catch up with 
the faster moving components. 

As light propagates in an optical fiber, some of it is 
scattered out of the fiber by Raman scattering leading to 
a loss of signal power. This is compensated for by the 
periodic insertion of short spans of lasing fiber. As the 
light passes through the lasing fiber, it is amplified by 
stimulated emission. It is impossible to avoid a certain 
amount of spontaneous emission in this process. This 
is amplified alongside the signal and provides a source 
of additive noise called amplified spontaneous emission 
(ASE) noise. 

Although we take full account of the effects of disper- 
sion, dispersion compensation and ASE noise in our anal- 
ysis, these are not the main focus this of work. A system 
with only these effects is a linear system with additive 
noise. This is precisely the-type of system understood 
by Shannon many years agoEJ. The additional aspects of 
optical communication are the non-linearities present in 
the propagation of light in an optical fiber. The most im- 
portant of these non-linearities is the optical Kerr non- 
linearity, by which the refractive index of the fiber de- 
pends upon the intensity of light in the fibenJ. 

The Kerr non-linearity causes scattering between dif- 
ferent frequency components of light in the fiber. Two 
incoming frequency components scatter into two outgo- 



ing frequency components. In WDM systems the effects 
of the Kerr non-linearity are grouped according to the ori- 
gin of the frequency components which scatter off one an- 
other. When the incoming and outgoing frequency com- 
ponents all lie in the same sub-band, the effect is known 
as self-phase modulation (SPM). It is this non-linearity 
that is responsible for the possibility of soliton propaga- 
tion within a sub-band. The combination of SPM and 
ASE noise leads to several mechanisms of signal distor- 
tion: fluctuations in the amplitude and arrival times of 
soliton pulses, known as Gordon-Haus jitterQ; non-linear 
amplification of the phase noisea; and, depending upon 
the sign of dispersion, a non-linear instability of ampli- 
tude fluctuations about a constant signal called the mod- 
ulation instabilityo. When the incoming frequency com- 
ponents are in different sub-bands and scatter back into 
their original sub-bands, the effect is called cross-phase 
modulation (XPM). This is the effect that we will concen- 
trate upon here. We will show explicitly that this leads 
to multiplicative noise- to a spatially and temporally ran- 
dom potential due to the signal in the other sub-bands. 
We will find a super-diffusive spreading of the pulse shape 
and a super-diffusive spreading of arrival times. Finally, 
the effect of the Kerr non-linearity when all of the incom- 
ing and scattered frequency components lie in different 
sub-bands is known as four wave mixing (FWM). From 
the point of view of the signal in any particular sub-band, 
it has an identical effect to additive amplifier noiseEj. 

A complete WDM optical communications system con- 
sists of a series of spans of standard, dispersive fiber in- 
terleaved with dispersion compensating fibers and loss 
compensating amplifiers. The precise ordering of these 
elements is known as the dispersion map. The perfor- 
mance of the system can depend quite sensitively upon 
this dispersion map. However, for concreteness, we will 
restrict ourselves to one particular map, shown in Figf. 

Our approach will be to calculate Feynman path inte- 
grals describing the propagation of light in a single sub- 
band, taking account of scattering from signals in the 
other sub-bands through XPM. The Green's functions for 
this propagation must be averaged over realizations of the 
signal in the other sub-bands. The Green's functions for 
propagation in the optical fiber turn out to be the Green's 
functions of a Schrodinger equation with spatially and 
temporally random potential. The main calculations in 
this paper are the evaluation of the single and two parti- 
cle Green's functions of this Schrodinger equation. This 
calculation is similar to calculations carried out in two 
other contexts: Turbulent flow or the growth of int erfaces 
may both be described by the KPZ equatior£Zrli§li3. This 
may be mapped to a diffusion equation with spatially and 
temporally random potential. Calculating the quadratic 
statistics for these processes and averaging over the tur- 
bulent flow or noise in the growth process is very similar 
to calculating the two particle Green's function in our 
problem. In particular, short-range temporal correlations 
allow the calculation to be reduced to an effective single 
body problem for which the Feynman path integral may 
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be readily calculated. We will find the same structure in 
our solution to the present problem. The second prob- 
lem to which our work bears similarity is the calculation 
of Green's function in the presence of a quenched ran- 
dom potentially. In this context, it is usual to truncate 
the diagrammatic series for the single and two particle 
Green's functions to the Born and Ladder series, respec- 
tively. Here we will find that this truncation is justified 
by the short range temporal correlations of the random 
potential. 

Our approach contrasts with other methods of includ- 
ing the effects of XPM in WDM systems], in which indi- 
vidual scattering events between soliton pulses in differ- 
ent sub-bands are considered. Our results ar,e very sim- 
ilar to those obtained for Gordon-Haus jittem, although 
of a different physical origin. 

Before going on to a detailed calculation, let us first 
state our results and the physics underlying them. As 
a signal pulse propagates along a particular sub-band it 
passes through regions of higher and lower intensity of the 
total signal in the other sub-bands. This causes the pulse 
to speed up and slow down due to the Kerr non-linearity. 
Since the signals in the other sub-bands are unknown to 
the users of any particular sub-band, this speeding up 
and slowing down of the pulse has a statistical uncer- 
tainty. The result is a diffusive spreading of the phase 
velocity and a corresponding super-diffusive spreading of 
the distribution of pulse arrival time; (5T 2 ) oc L 3 (to see 
this requires integrating the diffusively spreading phase 
velocity over the length of propagation). This is the main 
result of our calculation. The functional form is very sir 
ilar to that obtained previously for Gordon-Haus jitter 
although its physical origin is slightly different; it is due 
to XPM rather than a combination of SPM and ASE 
noise. In addition to the uncertainty in the pulse arrival 
time, the pulse shape itself is distorted due to the effect 
of XPM; since they travel at different speeds, the dif- 
ferent frequency components of the pulse sample slightly 
different portions of the signal in the other channels and 
arrive at slightly different times. Therefore, the pulse 
width itself will show a super-diffusive spreading due to 
the effects of XPM. The effect of XPM upon the pulse 
power is much weaker than its effect upon pulse timing. 
The processes leading to fluctuations in pulse power are 
of higher order in the (weak) non-linearity. The timing 
errors induced by XPM cause the pulse power to spread 
over a larger interval leading to the possibility of over- 
lap between adjacent pulses. This intersymbol interfer- 
ence allows timing errors to be converted into amplitude 
errors. Numerical simulation, however, shows that the 
dominant contribution to amplitude errors comes from 
SPMlIS, the effects of which are not studied here. 



II. THE MODEL 

The propagation of light in an optical fiber is described 
by Maxwell's equations with a dielectric constant depen- 
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FIG. 1: Schematic Diagram of Optical Fiber System 



dent upon the energy density of the electromagnetic field. 
The optical fibers used in communications systems sup- 
port only one transverse mode. Expanding the Maxwell's 
equation for the low frequency envelope of the electric 
field and taking the long wavelength limit, one obtains a 
non-linear Schr6dingeK-.equation with the roles of space 
and time interchanged!! 

id z E(z,t) = (/3d? -ia)E(z,t)--/\E(z,t)\ 2 E(z,t). (1) 

The term proportional to a on the right-hand side of 
this equation describes the loss of optical power from 
the channel and the term proportional to 7 encodes the 
effects of the Kerr non-linearity. Dividing the electric 
field into its components in each sub-band of the WDM 
system, the low frequency envelope of the electric field in 
the i th sub-band, Ei(z,t), is given by 

id z Ei(z,t) = ((3d? -ia)Ei(z,t) + V(z,t)Ei(z,t), 
V(z,t) = -27^^(2, t)\ 2 . (2) 

In writing down this expression we have neglected self- 
phase modulation, j\Ei(z,t)\ E i {z,t) 1 and four-wave 
mixing, ^Yj 3 ^i E *j( z ^) E k{ z ^) E i{ z ^)- Tnis neglect 
of self-phase modulation and four-wave mixing may be 
justified on two grounds. Firstly, there are regimes 
where numerical simulation confirms this to be a good 
approximation!!]. Studying the system in the absence of 
the effects of self-phase modulation and four-wave mix- 
ing allows us to deduce the effects of XPM. Our aim is to 
determine what constraints are imposed upon communi- 
cation by these effects. In a real system, other sources of 
signal distortion will impose additional constraints, but 
these are not of direct concern in this paper. Equation^) 
corresponds to a non-linear Schrodinger equation with a 
spatially and temporally random potential due to non- 
linear interaction with the signals in the other sub-bands. 
The effects of loss in this system may be accounted for by 
using rescaled electrical fields that maintain their ampli- 
tude during propagation; E e ff t i(z,t) = e az Ei(z,t). The 
propagation of these rescaled fields is described by Eq. (j|) 
with an exponential rescaling of the strength of non- 
linearity, j e ff = j(l — e~ aL )/aL, where L is the length of 
a single span of fiber, and without the loss term propor- 
tional to a. Henceforth, we will assume the use of these 
rescaled fields and, drop the subscript eff for brevity. 

We follow Refills and model the stochastic potential by 
a Gaussian distribution with short-ranged correlations in 
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space and time; the signals in the separate sub-bands are 
statistically independent and have short-ranged temporal 
correlations on timescales of order 1/bandwidth. Rel- 
ative dispersion of signals in different sub-bands leads 
to short-ranged spatial correlations, which may be ap- 
proximated by a delta-function if the dispersion is suf- 
ficiently large. The strength of the potential, rj = 
fdz(SV(z,t)SV(Q,t)) = [( 1 P) 2 /{2(3A5W)]\n[n c /2], is 
obtained by summing the contributions from n c channels 
with separation 6W, bandwidth A and signal power, P. 
The logarithmic dependence upon the number of chan- 
nels is due to the suppression-of the effects of widely sep- 
arated channels by dispersior£9. In the following section, 
we will often write (5V(z,t)SV(0, 0)) = 2nr]/A5(t)5(z), 
with an explicit delta-function in time. This simplifies 
our calculations and elucidates which contributions are 
most important for propagation in a potential with short 
temporal correlations. It is important to realize, however, 
that the potential and signal are both band-limited, so 
that the correlations are not strictly delta-functional. We 
account for this by regularizing the delta-function of zero 
argument to S(t = 0) = A/2n. 

The system considered here consists of 10 x [80km 
standard fiber followed by 16km dispersion compensat- 
ing (d.c.) fiber and a loss compensating amplifier]. The 
input signal has power, bandwidth and channel sepa- 
ration P = 5mW, A = 10GHz and SW = 15GHz, 
respectively. The standard/dc fiber loss, non-linearity 
and dispersion parameters are a — 0. 048/0. 115km -1 , 
7 = 1.2/5.1W -1 km _1 and = 11.0/ - 66.9ps 2 km _1 . 
Comparing jP/a for the standard and d.c. fibers gives a 
measure of their relative strengths of non-linearity. The 
power entering the d.c fiber is a factor of ~ 50 lower 
than that entering the standard fiber. This, combined 
with the greater loss in the d.c. fiber, leads to an effec- 
tive strength of non-linearity ^150 times greater in the 
standard fiber. We ignore non-linearity in the d.c fiber 
in our analytical work. It may be included in precisely 
the same way as in the standard fiber and does not lead 
to any qualitative change in our results. Three impor- 
tant lengthscales may be determined from these system 
parameters; the total fiber length, L to t = 10 x 80km, 
the non-linear length, r\~ x = 223km, and the dispersion 
length, L D = (/3A 2 ) -1 = 747km. The response of the 
system is completely determined by the ratios of these 
three lengthscales. 

III. INPUT OUTPUT STATISTICS 

The output electric field, E out (t), after propagating 
along a single span of standard fiber, may be expressed 
in terms of the input electric field E in (t) as 

E out {t) = J dt'g(t,t')E m (t')+n(t), (3) 

where n(t) is additive Gaussian white noise, with vari- 
ance N, representing ASE noise. Q(t,t') is the Green's 



function of Eq.(g) for propagation from z = at time t 
to z = L at time t' , with a fixed realization of V{z,t). 
More generally, we will use the notation Q(z%, z%\ ti, *2) 
to indicate the Green's function for propagation from 
t = t\, z — z\ to t = t2, z = Z2- Using this notation, 
Q{t\,t2) = Q(0, L; ti, t2). Suppression of one of the spa- 
tial indices implies propagation from z = 0; Q(z; t\, t 2 ) = 
Q{0,z;t\,t 2 ). Similarly, Q{z,t) = 0(0, z;0,t). It is par- 
ticularly important to note that Q(t\, t 2 ) ^ Q{t\ — t 2 ) as 
would be the case for a static potential. 

Dispersion compensation is included by replacing 
Q(t,t') in Eq.(||) with a dispersion compensated Green's 
function, Qd{t,t'). This is obtained by convolving the 
Green's function for the standard fiber, Q(t,t'), with 
the Green's function for the dispersion compensating 
fiber, Q comp (t,t'); G d {t,t') = J dt" Q comv (t,t")Q(t" ,t'). 
Gcom P (t,t'), has the opposite sign of (3 to Q{t,t') and de- 
scribes propagation with no potential, V . 

Propagation along a series of interleaved standard and 
d.c spans is modeled by convolving a string of disper- 
sion compensated Green's functions. This convolution 
is independent of the ordering of the Green's functions 
and the effects of XPM are, therefore, independent of 
the dispersion map. Real systems are sensitively depen- 
dent upon the choice of map due to the effects of SPM, 
which is ignored here. Finally, due to the unitarity of 
Qd(t,t'), adding ASE noise at each amplifier is equiva- 
lent to adding Gaussian white noise at the receiver. An- 
alytical calculations may, therefore, be performed using 
a simplified dispersion map consisting of N s standard 
spans followed by N s dispersion compensating spans and 
adding Gaussian white noise at the receiver. 

We now use this simple model of the system to calcu- 
late the effects of XPM upon the input/output statistics. 
The quadratic statistics are given by 

(E* n (t)E m (t')) = P6(t-f), 
(E* n (t)E out (t>)) = P(G d (t,t')), 
(E: ut (t)E out (t<)) = (P + N)S(t-t'), (4) 

and the quartic statistics are given by 

(S\E ln (t)\ 2 S\E ln (t')\ 2 ) = AP 2 5(t-t')/2ir, 
(S\E m (t)\ 2 S\E out (t')\ 2 ) = P 2 g I d I (t-t>), 
(5\E out (t)\ 2 6\E out (t')\ 2 ) = A{P + N) 2 5{t-t')/2^(b) 

The angular brackets indicate averages over the signal 
in the sub-band of interest, the signals in the other sub- 
bands and the amplifier noise. {Gd(t, t')) is the average of 
the dispersion compensated, single particle Green's func- 
tion over the potential (or realizations of the signal in 
the other sub-bands) and Q^(t - t') = (\Gd(t, t')\ 2 ) is 
the average of the two-particle Green's function over the 
random potential. The input signal is assumed to have 
a Gaussian distribution with power P and the ampli- 
fier noise to have a Gaussian distribution with power N . 
We have used the unitarity of the single point Green's 
function for a particular realization of the potential in 
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deriving Eqs.(§) and (||); Q{t,t) = 1. The delta-function, 
6(t—t' = 0), has been regularized as discussed in Sec. II, to 
allow for the fact that the sub-band has a bandwidth A; 
S(t — t' = 0) = A/27T. The averaged Green's functions, 
(Gd{t,t')) and Gd (t — t') are invariant under temporal 
translations and are, therefore, functions oft — If. 

The statistics embodied by Eqs.(H) and (0) are non- 
Gaussian, since G^ 1 (t - f) + \{\G2tJ))f in general. 
The corrections to Gaussian statistics are given by the 
vertex corrections to the two-particle Green's function. 
Alternatively, one may say that linear propagation in the 
presence of a random potential leads to non-linear prop- 
agation on average. Calculation of the quadratic and 
quartic statistics in terms of Green's functions allows a 
simple characterization of this non-linearity. In the fol- 
lowing section, we will give an explicit calculation of the 
averaged Green's functions and follow this in the subse- 
quent sections by a discussion of their consequences for 
optical communication. The quadratic statistics will be 
useful for our discussion of coherent communication and 
the quartic statistics will be useful for our discussion of 
incoherent communication. 



IV. CALCULATIONS 



The average Green's functions (G) and Q 11 may 
be calculated ■— from their Feynman path integral 
representationsli-l If we assume that the stochastic po- 
tential is delta- function correlated in time, these path 
integrals may be evaluated in closed form. These results 
correspond to re-summing the Born and ladder seriesE-3 
in the diagrammatic expansions of (G) and Q 11 , respec- 
tively. These series may be used to approximate the ef- 
fects of a non-delta-function correlated potential. The 
diagrams that arc neglected in this approximation are 
small for a potential with short ranged temporal correla- 
tions; the use of a delta-function correlated potential is 
a calculational tool that emphasizes the range of validity 
of this approximation. In the following , we will calcu- 
late equations of motion for (Q) and Q 11 from their path 
integrals, under the assumption that V is delta-function 
correlated in time. We will show how the solutions of 
these equations correspond to the re-summed perturba- 
tive, ladder and Born series. Finally, the average Green's 
functions (G) and G 11 will be calculated allowing for the 
band limit of the sub-band and potential due to the other 
sub-bands. 



A. Single particle Green's function 

Recalling that the roles of space and time have been in- 
terchanged here compared with their more familiar roles 
in quantum mechanics, the single particle Green's func- 



tion is written in the form 

,t(L)=t 

(G(L,t)) = ( Dt(z)exp 

Jt(0)=0 

-i { dzV(z,t(z)) 
Jo 



t{L)=t 



t(0)=0 



Dt(z) exp 



dz±(d,t(z)f 

Z P 



dz^-(d z t(z)) 2 
o Z P 



L r L 



2 Jo Jo 



dzdz'(V{z, t(z))V(z', t(z')))(§) 



where the angular brackets denote averages over realiza- 
tions of V. We have carried out the average over V in 
order to obtain the second line, assuming that V is a 
Gaussian random variable. _ 

Taking the derivative of Eq. (^) with respect to Ltd, we 
obtain the following equation of motion for (Q): 



-i^d\-d h \ (c7(I ./)) 
= 6(L)6(t) 



dt'dz{V(L,t)V(z,t')) 



x(G(L-z,t-t'))(G{z,t')). 



(7) 



This equation is local in time for a stochastic 
potential that is delta-function correlated in time; 
{V(z,t)V (0,0)} = (2wr]/A)S(t)S(z). The solution of this 
local equation of motion is (G(L,t)) = e~ ,lL / 2 Qo(L,t), 
where Go(L,t) is the Green's function of the operator 
(— ipdf/2 — 8l), i.e. the Green's function for propaga- 
tion^ the absence of V. This result was first obtained in 
RefM Equation® may also be obtained diagrammati- 
cally by re-summing the Born series; it is nothing but the 
Dyson's eauation for (Q(L,t)) obtained by re-summing 
this seriesEj. The temporal delta-function correlation of 
V implies that the propagators of V must not cross, thus 
restricting the contributing diagrams to the Born series. 
Figure 2. shows a symbolic representation of Eq.® after 
multiplying on the left by Go- The dashed line indicates 
propagators of V, the solid line indicates £o and the thick 
solid line indicates (G). This diagram is to be interpreted 
in real space or momentum space in the usual way. The 
momentum/frequency space solution is 



(G( P ,uj)) =i(p- (3uj 2 + ir]/2y 



(8) 



where p is the wavevector conjugate to z and uj is the 
angular frequency conjugate to t. For later calculations, 
it is useful to have (G) in a mixed, position/frequency 
notation, where it is given by 



(G{z, u)) = exp [- (77/2 - ipoj 2 ) z] 



(9) 



Dispersion compensation is achieved by temporally con- 
volving G(L,t) for the standard fiber with that of the 
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FIG. 2: Born series for (Q). 

compensating fiber (or alternatively multiplying together 
the position/frequency space Green's functions). Ig- 
noring non-linearity in the compensating fiber, we find 
Gcom P {L,t) = C/q(L,£). The dispersion compensated 
Green's function is then given by 

(g d )(L,t)=e-^ 2 6(t)- (10) 

On average, therefore, the effects of XPM non-linearity 
separate from those of dispersion. Moreover, the effect 
of XPM is to dephase the propagating electric field. 

B. Calculation of Q u 

The calculation of Q 11 from its path integral proceeds 
in a very similar manner to the calculation of (Q). Af- 
ter writing down the path integral representation of Q 11 , 



we perform the average over V and deduce an equation 
of motion for Q H by differentiating with respect to L. 
This differential equation is local in time when the cor- 
relation function of V is a temporal delta-function. The 
equation of motion is most easily solved in the frequency 
domain, where it is recognizable as the re-summed lad- 
der series. As before, the set of contributing diagrams is 
restricted due to the requirement that the propagators of 
V may not cross. We will first of all derive an expression 
for Q IJ (u}, q) in the absence of dispersion compensation 
from its path integral representation. We will verify that 
this result is identical to that obtained by summing the 
ladder series of diagrams. In fact, closed expressions for 
<7 /J (-L,i) are most readily obtained by re-summing this 
series directly in the length/frequency domain. We ex- 
plicitly carry out this re-summation to obtain expressions 
for Q 11 (L, uj) and its dispersion compensated counterpart 

We will first derive a path integral and equation 
of motion for the Green's function, Q 11 (L;ti,t 2 ) = 
{G(L, Ti)Q*(L, T2)). In the end we will be interested in 
g H {L,t) = £ /7 (L;ti = t,T 2 = t). The path integral 
representation of Q 11 is 



g u (L;n,T 2 ) = {g{L lTl )g*{L,T 2 )) 

l>Ti (L)=ti ,t 2 (L)=t 2 



Dti(z)Dt 2 {z) exp 



ri(0)=0,r 2 (0)=0 

1 

2 

t 1 (L)=t 1 ,t 2 (L)=t 2 



dz 



(d z n) 2 , {d z r 2 f 

2(3 % 2(3 



x exp 



dz (-W (z,n(z)) +iV (z,t 2 (z))) 



T!(0)=0,r 2 (0)=0 
1 

x exp 
x exp 
x exp 



Dt 1 {z)Dt 2 (z) exp 



dz 



(d z n) 2 .(8 z t 2 ) 
2/3 1 2/3 



2 \ 1 



- J dzdz'iViz.r^Viz'.niz 1 ))) 
~ [ dzdz'(V (z,t 2 (z))V (z',t 2 (z'))) 

dzdz' \V (z,Ti(z))y (z',T 2 (z'))> 
I 



(11) 



The average over V has been carried out in passing from 
the penultimate to final line of Eq. ( |i~i] ) . The equation of 
motion for g 11 is derived by differentiating both sides of 
this equation with respect to L. The resulting differen- 
tial equation is local in time when V has delta-function 
correlations in time. It is slightly easier to make the as- 
sumption of temporally delta-function correlations for V 
in Eq.(|ll|) before differentiating. The path integral then 
reduces to 



x exp 
2tt 



dz 



(d z n) 
2(3 



. (d z r 2 ) 
' 2/3 



V-^5( Tl {z) -r 2 {z)) 



+ 1] 
(12) 



g u {L- Tl ,T 2 ) = 



Tl(L)=Tl,T 2 (L)=T 2 



Tl(0)=0,T 2 (0)=0 



D Ti {z)Dt 2 {z) 



This result is very similar to that obtained in the analysis 
of the diffusion equation with a spatially and temporally 
random potentialO. In that case, the path integral for 
the square of the wavefunction reduced to that of two 
interacting random walks. Here the path integral for g 11 
reduces to that of two interacting Schrodinger particles, 
or waves, propagating in opposite directions. Differenti- 
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FIG. 3: Ladder series for Q 



ating Eq.(jT^) with respect to L, we obtain the following 
local equation of motion for Q 11 . 

= <5(£)^(r 1 ) ( 5(r 2 )+7 ? J ( 5(T 1 -r 2 )e // (L;r 1 ,T 2 )(13) 

This equation is most easily solved by first changing co- 
ordinates to the mean time, t = (t± + t 2 )/2, and time 
difference, T = (n — t 2 )/2. After Fourier transforming 
in space and time we find 

-i(q- + w?£(0)^ G n {q; u, fi) 

= l + ^/fV(^,Q), (14) 

where a; and 17 are the angular frequencies conjugate to t 
and T, respectively and q is the wave-vector conjugate 
to L. Equation (O) is equivalent to the ladder series 
shown in Fig. 2. In order to make this correspondence 
clear, the effects of non-linearity have been divided into 
two parts; the term proportional to r\ on the left hand 
side of Eq. ( |l4| ) takes account of self-energy correction to 
the single particle Green's function and the term propor- 
tional to r\ on the right hand side takes account of vertex 
corrections. Following the usual convention, we denote 
the solution of Eq.fJlJ) in the absence of vertex correc- 
tions by n(<7; u, f2) = —i(q — 0u>Q,/2 + irj) . One may 
check by direct substitution of (Q(q, lj)) from Eq.(||) that 
U(q;uj,n) = f dp/2Tr{g(q+p,uj + n))(g*( P ,n)). Finally, 
using the notation 11(3, a; ) = / dfl/27rll(q; to, f2), the so- 
lution of Eq. (O) is 



g u {q^) 



n(«,«) 



1-^nfoa,) 



E(g,w) [l 



2lir\ 



n(g,w) 



(15) 



appropriate diagrammatic series is shown in Fig. 4. The 
result of summing this series is 

= IU(L,u) 

+—^- J dxU{x,Lu)U d (L - x,w) 

+ I dX J d y^ x ^) n ^^)^d{L-y^) 

+... (16) 
n(z, uj) is given by 

n(z,u) = J^(g(z,u/ + w/2))(g*(z,<J-w/2)) 

= e-/^- (17) 

and 11^(2,0;) is the dispersion compensated analogue of 
n(z,£j), given by 

ITi(z, W ) = J ^-(g{z,u' + uj/2))g* {L,uj' + uj/2) 

x (0*(z, a/ - u/2))g (L, - w/2) 
= / ^_ e ip u ' u (»-L) M (18) 

J 27T 

The calculation of g 11 proceeds by substituting LT^z, u>) 
and II(z, w) into Eq.(|lq), carrying out the frequency inte- 
grals in each term, followed by the integrals over length, 
and finally summing the contributions from each term. 
In carrying out this procedure, one must take careful 
account of the support for the various frequency inte- 
grals. The band limitation of the sub-band implies that 
the angular frequency carried by each internal Green's 
function must lie between — ttA and +7rA. Carrying 
out the frequency integrals with this restriction leads to 
very complicated expressions. An alternative approach, 
which yields much more manageable expressions, is to 
take into account the fact that the stochastic potential 
is not strictly delta-function correlated in time, but is 
band-limited to lie between — ttA and +ttA. We use the 
fact that V has short range correlations in time to justify 
neglecting terms other than those of the Born and ladder 
series in our diagrammatic calculation, but allow for the 
band limitation of V in calculations of internal frequency 
integrals. 

Let us illustrate this by explicitly calculating the sec- 
ond order term in Eq.(|16[). This is given by 



which corresponds to the ladder diagram shown in Fig.3. f dx [ dyU(x, (j)U(y, uj)U d (L - y, uj) 

Jo Jo 



The modification to these results due to dispersion 
compensation is not quite as simple as the modification 
to (g). Although it is possible to derive a path integral 
expression for akin to Eq. (|l3|) , it is easiest to derive 
the dispersion compensated g 1 Rising diagrammatic per- 
turbation theory in the length/frequency domain. The 



, f x , f dh)\ duj2 duj 3 

dx dy 

o J y J 2tt 2tt 2tt 

x exp [if3ui(uj3 — L02)x + if3uj(uj2 — wi)y] 



A 



2W Jo ,/n X Y 
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standard fiber 



d.c. fiber 





n(co, x) n(co,y-x) : n d (co,L-y) 



FIG. 4: Position-space diagrammatic series for Q^ 1 (L,i 



2ir 2 I 



(19) 



X=ki3ujAL 



where SifX] = J Q dusmu/u. We have carried out 
the frequency integrals in passing between the second, 
third and fourth lines. The integration variables have 
been changed from w\, u>2 and 0*3 to u>i, L03 — u>2 and 
LU2 — lo\ ■ Integrations have been carried out over the do- 
mains —A/2 < wi,W3 — 0^2,^2 — < A/2, the first 
corresponding to the band-limitation of the signal in the 
sub-band and the second and third to band-limitation of 
the stochastic potential leading to a restriction on the 
transfer of frequency at each vertex. 

Calculating each term in the series in this way and 
re-summing, we obtain our final result 



- V L 1 - 



Si\X] 
X 



X=tt(3ujAL 



(20) 

for dispersion compensated propagation. The Green's 
function for uncompensated propagation is calculated 



similarly replacing 11^ with II in Eq.(|16[). The result is 



g"(L, u ) = g 1 d 1 (L,u) 



sin A 



X 



(21) 



X=k(3ujAL 



The prefactor of sin A/A describes the dispersive spread- 
ing of a pulse propagating in the absence of non-linearity. 
The distortion of the pulse due to XPM is described by 



Q^(L,uj). There are two important limits of Eq.(|20|) 



At small frequencies, Si[A]/A — > 1 - X 2 /18. As a 
consequence, ^ j (£,cj) tends to a Gaussian of width 
v/p) = 27rA e// = (3A/n)^Ll/ V Ll t . The interpre- 
tation of this width as an effective bandwidth will be 
discussed below. At large frequencies, Q¥{L,uj) tends to 
a constant. These limits imply that Q d (L,t) has Gaus- 
sian tails at large times and approaches a delta-function 



at short times. We will discuss the consequences of these 
results for communication in the following sections. 



V. AMPLITUDE AND TIMING JITTER 

The results of sections III and IV give the input /output 
statistics of the WDM optical fiber in terms of its single 
and two-particle Green's functions. In this and the fol- 
lowing section, we shall use these results to discuss the 
effectct of XPM upon communication. 

Firstly, we discus two operational measures of the ef- 
ficiency of an optical communications system. Our dis- 
cussion will be quite brief, since these considerations are 
not central to the main message of this paper. Neverthe- 
less, it is important to demonstrate that analytical results 
such as Eqs. © and @ can be used to enhance under- 
standing from this operational perspective. A common 
way of characterizing signal distortion in optical com- 
munications is by amplitude and timing jitter. These 
measures of system performance are specific to digital 
communication. A digital message is made up of a string 
of pulses or marks representing Is, and spaces represent- 
ing zeros. As a pulse propagates along a noisy channel, 
it speeds up and slows down in a random way leading to 
a distribution of arrival times(defined as the centroid in 
time of the power distribution at the output) at the re- 
ceiver known as timing jitter. Similarly, the total power 
of the pulse fluctuates leading to a distribution of am- 
plitudes of the received pulse known as amplitude jitter. 
The distributions of amplitude and timing fluctuations 
need not be independently. In fact, the terms amplitude 
and timing jitter are often reserved for the variance in 
amplitude and arrival times of a pulse. This is the sense 
in which we shall uses these terms from now on. In ad- 
dition to jitter, the shape of the individual pulses may 
be distorted during propagation. We do not discus these 
distortions here. 

The effects of SPM and ASE upon the amplitude and 
timing jitter of soliton, or soliton-like, pulses has n been 
considered in a number of works, starting with Ref.Q. All 
of these works expand in small fluctuations- induced by 
the additive ASE noise- about a soliton solution. Such 
analyses lead to the conclusion that amplitude jitter is 
not very significant and that timing jitter leads to an L 3 
dependence of the variance of arrival times. This latter 
result iti known as Gordon-Haus jitter affects discovery 
in Ref.Q. A recent work of Falkovich et a&3 has used a 
saddle point approximation within a Martin-Siggia-Rose 
field theory for this propagation to obtain a full joint 
distribution of amplitude and timing errors. This work 
shows in a particularly transparent way, the connection 
between the method of optimal fluctuations in field the- 
ory and large deviations in statistics. Here we use Green's 
function techniques to consider the contribution of XPM 
to jitter. Despite its different physical origin, the result- 
ing jitter has a broadly similar dependence upon system 
parameters to Gordon-Haus jitter. 
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In principle, the output for an arbitrary input pulse 
shape may be calculated, accounting for the effects of 
XPM, in terms of Green's functions for propagation along 
the fiber. The amplitude and timing jitter may then be 
obtained by taking appropriate averages. This calcula- 
tion may in practice prove rather cumbersome. How- 
ever, the functional dependence of the jitter upon sys- 
tem parameters is expected to be broadly independent 
of the input pulse shape (the magnitude of the jitter will 
show a weak dependence upon the pulse shape). For 
simplicity, therefore, we present approximate expressions 
for amplitude and timing jitter due to XPM, by con- 
sidering the output for constant and delta-function in- 
puts, respectively. The simple expressions obtained in 
this way demonstrate the functional dependence of the 
jitter upon the system parameters. In particular, we find 
that the variance of arrival times is proportional to L , 
whereas amplitude jitter is largely independent of XPM 
. These functional dependences are the same as those 
due to Gordon-Haus jitter although the physical origin is 
somewhat differentia. 

Timing jitter may be found by considering the re- 
sponse to an input delta-function in time. A real sig- 
nal pulse is band-limited, but it may be thought of as 
being made up of many delta functions. The use of 
a delta-function input allows us to calculate the func- 
tional form of the timing jitter in an analytically tractable 
way. The dependence of the resulting jitter upon sys- 
tem parameters is expected to be broadly similar for 
arbitrary input pulse shapes. Without XPM, the out- 
put electric field resulting from a delta function input is 
E ou t(t) ~ Gd(t) oc S(t) and the modulus of the electric 
field is given by \E out (t)\ 2 = \Qd{t)\ 2 oc 5(t). Introduc- 
ing XPM broadens this response. Different frequencies 
making up the pulse see a slightly different stochastic 
potential and arrive at slightly different times. Since the 
signals in the neighboring sub-bands are unknown, this 
broadening leads to uncertainty in the arrival time. Av- 
eraging the response to a delta function input over real- 
izations of the signals in the other sub-bands, therefore, 
we find that Q 1 / {t) = {\Gd(t)\ 2 ) gives a measure of the 
distribution of arrival times. Using this interpretation 
and the results of the preceding section on the calcula- 
tion of G 1 /, we find (t 2 ) = {2nA eff y 2 oc L^ ot rj; XPM 
causes a super-diffusive spreading of arrival times. This 
functional dependence of the timing jitter upon the fiber 
length is identical to that found by Gordon_and Haus for 
the combined effects of ASE noise and SPMO. Notice that 
this contribution to timing jitter arises entirely from ver- 
tex corrections to the 2-particle Green's function, Q^ 1 (t); 
if we ignore vertex corrections by making the approxima- 
tion Q^(t) = \ (G d (t)}\ 2 and substituting Eq. @, we find 
no contribution to timing jitter. 

Using Eq.(|3|) to calculate the statistics of the output 
for inputs held constant at Ei n (t) — and Ei n (t) = A 
gives an estimate of the amplitude statistics of a I and 
in the output bit stream!! For a zero input, the output 
has the same Gaussian distribution as the noise. For an 



input of constant amplitude A, 

(Eoutit))! = A{g d (u = Q)), 
(\E out (t)\ 2 ) 1 = (A 2 +N)S(0), 
(S\E out \ 2 (t)S\Eout\ 2 (t))i = A i [(\Q d (uj = 0)\ i )-5(0) 2 } 

+N(2A 2 + N)6{0) 2 
w N(2A 2 + N)6(0) 2 , (22) 

where the delta-function of zero argument is to be un- 
derstood as A/27T = (5(0) for our band- limited sub-band. 
The evaluation of these results is outlined in the Ap- 
pendix. The output for a 1 input is not Gaussian: a 
Gaussian fit to (E out (t))\ and (\E out {t)\ 2 )i does not re- 
produce (5\E out \ 2 (t)6\E out \ 2 (t))i. These results also pre- 
dict that the contribution of XPM to amplitude jitter 
is small. Equation ( p2| ) depends upon the evaluation of 
(|£?(w = 0)| 4 }, which we have not carried out here. A 
full evaluation of this 4-particle Green's function would 
be required to give a full expression for amplitude jitter. 
If, however, we approximate (\Q(u> = 0)| 4 ) « (\Q(u) = 
0)| 2 ) 2 = <5(0) 2 then amplitude jitter shows no dependence 
upon XPM. Corrections to this result arise from higher 
order vertex corrections to the 4-particle Green's func- 
tion and are small for weak non-linearity. In particular, 
these vertex corrections are of higher order than those 
that give rise to timing jitter. The predominant effect of 
XPM, therefore, is to induce timing errors. This conclu- 
sion is in accord with other works which have attempted 
to model the effects of XPM as a modification to the 
Gordon-Haus results for soliton jittem. In these works, 
individual scattering events between soliton-like pulses in 
different sub-bands of the WDM system are considered 
and the amplitude jitter due to scattering ignored. 

VI. INFORMATION THEORY 

In this section, we will discuss how the input/output 
statistics calculated in sections III and IV may be used 
to construct general arguments about the maximum rate 
of communication via a WDM system. This is still an 
open problem. The only systems for which exact answers 
about information capacity may be obtained are those 
with Gaussian input/output statistics. We will use the 
results for Gaussian systems to obtain lower bounds upon 
and approximations to the capacity of the present, non- 
Gaussian system. 

The ideas of information theory propounded by Shan- 
non are very much akin to the ideas of statistical me- 
chanics. For a signal made up of a sequence of let- 
ters (or field values), chosen from some alphabet (or 
range of values) one may define a probability distribu- 
tion, p(x), for the probability that a particular letter 
takes the value x. The amount of information per let- 
ter in such a sequence is given by the entropy of this 
distribution; H{X) = H[p(x)} = — '}2 x p(x)\ogp(x). H 
the logarithm is taken in base 2, this gives the amount 
of information per letter in bits. One way to understand 
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this is to think of compressing a signal. Compression 
involves removing redundancy in the signal (or alterna- 
tively correlations between letters) and so increasing the 
amount of information per letter. In this way, a signal 
with higher entropy will carry a higher density of infor- 
mation. 

When one attempts to communicate via a channel, 
each letter, x, of the input signal will be received as 
some letter, y, at the output. Allowing for the uncer- 
tainty in the transmission process, the output for an in- 
put, x, will be given by some conditional distribution 
p(y\x). For an information theorist, this defines the chan- 
nel of communication. For each input letter, the out- 
put may take a range of values. The entropy of this 
additional spreading at the output is given by the con- 
ditional entropy; H(Y\X) = J2 x p( x ) H ( Y \ X = x ) = 

- J2 X p( x ) E„ p(v\ x ) k>gp(y\x)] . 

When a sequence of letters is transmitted, the amount 
of information in common between the input and output, 
otherwise know as their mutual information is given by 
the total entropy of the output signal minus the addi- 
tional entropy due to signal distortion; 



spectral efficiency, is given by integrating over frequency; 



I(X;Y) = H(Y) - H(Y\X) 



(23) 



When using a particular channel, the conditional distri- 
bution p{y\X) is fixed and the remaining freedom is in 
the choice of input distribution, p(x). The idea is that 
if certain letters are more distorted than others, it pays 
one to use them less frequently, even though this means 
a reduction in the amount of information that the input 
signal can carry. The maximum transmitted information 
density that can be achieved is given by functionally op- 
timizing the mutual information over the choice of input 
signal distribution; 



C = max I(X;Y). 

p(x) 



(24) 



Shannon formulated these ideas rigorously and showed 
that this is not only an upper bound upon the rate of 
transmission of information, but that it is also an achiev- 
able bound. A summary of.±hese basics of information 
theory may be found in Ref.0. The functional optimiza- 
tion in Eq.(j2^) is carried out under the various system 
constraints, such as the average signal power. In fact, 
the only case for which there exists an explicit analyti- 
cal solution is when the joint distributions of input and 
output signal are Gaussian. In this case, the capacity is 
given by 



C* = loe 



(x*x)(y*y) 



(x*x) (x*y) 

(y*x) (y*y) 



(25) 



When the signal is band-limited, the corresponding ex- 
pressions for the capacity per unit bandwidth, or the 



2ttA 



dui log 



(y>-uj) {yZv- u ) 



(26) 



The factor of 1/2tt arises because u> is the angular fre- 
quency conjugate to t. 

We are now in a position to use these results from 
information theory, together with our knowledge of the 
quadratic and quartic statistics of the WDM system to 
obtain a lower bounds and estimates of the capacity of 
the WDM optical communications system. A rigorous 
lower bound on the capacity of a channel defined by a 
conditional distribution, p(y\x), may be found as follows: 
The quadratic input/output statistics are calculated for 
a Gaussian input distribution. A channel whose statistics 
are Gaussian and whose quadratic statistics correspond 
to those of the actual channel, will have capacity given by 
Eq. ( ^6|) . This gives a lower bound on the capjacity of the 
actual channel. A proof of this is given in Ref.li3. The ba- 
sic idea behind it is that deviations from Gaussian statis- 
tics correspond to extra correlations in the noise that the 
channel adds to the signal. These extra correlations al- 
low one to make progress in determining which part of 
the received signal is noise and which part is signal. Ig- 
noring these correlations- as one does in the Gaussian 
approximation- leads one to underestimate the extent to 
which noise may be removed from the signal and so to an 
underestimate of the system capacity. Notice that since a 
Gaussian distribution maximizes the entropy for a fixed 
variance, any deviation from Gaussianity will increase 
correlations. 

We are interested in two ways of using the WDM op- 
tical fiber system. In the first case, both the amplitude 
and phase of the electrical signal are used to communi- 
cate. We call this coherent communication. In the sec- 
ond situation, only the amplitude of the electrical field is 
used to communicate. This latter case is more represen- 
tative of current optical communication systems. Coher- 
ent coiainunication was considered by Mitra and Stark 
in Ref.Ej. The quadratic statistics given by Eq.([|) are 
combined with Eq.(^6|) to give the following expression 
for the spectral efficiency in terms of the single particle 
Green's function of the channel: 



c > 



2ttA 



log [1 + P ef f(u)/N ef f(u)] 



(27) 



with the effective signal and noise powers given by 



p e// H = p|(^M)| 2 

N eff (u,) = P (l - KOiM)! 2 ) + N, (28) 

respectively. This result is a modification of the famous 
ShannonEJ result for the capacity of a linear channel with 
additive white noise, allowing for the conversion of sig- 
nal power to noise power by scattering from signals in 
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the other sub-bands. The single-particle Green's function 
calculated in section 4, (Q d (L,uj) oc exp[— const, x P 2 L] 
is substituted into Eq.(|8|) to obtain the functional de- 
pendence of the capacity upon system parameters. The 
result is sketched in Fig. 5. There are two main features 
worthy of comment. The first is the peak in capacity 
for a particular power; increasing the signal power rela- 
tive to the noise power is productive up to a point, but 
when the signal power is very large, non-linear interac- 
tions between signals in the different sub-bands become 
the dominant mechanism of signal distortion. The ap- 
pearance of an optimal power is a common feature of the 
effects of noniJinearity. It occurs in the Gordon-Haus jit- 
ter of solitojuaa and in the non-lineaj-/Gordon-Mollenauer 
phase noiseO(see for example Ref.E3). The second feature 
worthy of comment is the dependence of capacity upon 
length at large distances. This is important, since it is for 
very long distance communication that propagation non- 
linearities have the most important effect. The Gaussian 
approximation to coherent communication gives a lower 
bound on capacity that decays exponentially with length 
at large distances. We do not expect this exponential 
decay to correctly reproduce the dependence of the sys- 
tem capacity, since it is in this limit that non-Gaussian 
corrections to the input/output statistics are likely to be 
most important. We will see below how taking account 
of these non-Gaussian corrections is likely to change the 
dependence of capacity at large distances. 

Unfortunately, there exist no rigorous bounds on the 
capacity for incoherent communication in which the 
phase of the electric field is ignored. Nevertheless, ex- 
pressionssimilar to Eq.(28) can provide useful estimates 
of the incoherent capacity. One such estimate is obtained 
by fitting an approximate Gaussian distribution of in- 
tensity about its mean values (given by Eq.(3)) to the 
quadratic statistics of intensity given by Eq. (| ) . The ca- 
pacity of the channel with this Gaussian distribution is 
an approximation to the capacity of the actual channel. 
It is given by 




2Tr\g I d 1 (u;)\ 2 /A) + 2PN + N 2 " 



The pre-factor of 1/2 arises because the intensity is 
a real field. In the limit of strong non-linearity and 
large signal to noise ratio, Eq. (pSj) reduces to C w 
v /r nj2A eff U 3/2 [P 2 /(P + N) 2 ]. This has an appealing 
interpretation. The effects of non-linearity are contained 
entirely within a prefactor and provide an effective band- 
width A e ff = l oAyjL 2 D lr\L\ ot . Our analysis of jitter 
showed that XPM mainly introduces timing errors and 
not amplitude error. The effective bandwidth enters be- 
cause signal pulses must be separated by more than the 
timing error in order to be distinguishable. The infor- 
mation carried by each pulse is determined by the ad- 
ditive noise power. Notice that in the limit of strong 
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FIG. 5: Spectral Efficiency vs Input Power: Plots of i. 
Eq.(|2^)-solid black line; ii. Shannon formula, c = log[l + P/N] 
-short-dashed, red line; iii.Eq.(29)- long-dashed, green line. 



non-linearity, the coherent capacity has an exponential 
decay with length, whereas Eq.(p9|) has a power-law de- 
cay. The actual coherent capacity is greater than the 
incoherent capacity, since a greater number of degrees of 
freedom are used in coherent communication. The be- 
havior of Eqs. (|2^) and ( p9|) at large powers point to lim- 
itations of the Gaussian approximation. A comparison of 
the coherent and incoherent capacity is shown in Fig. 5. 
The incoherent capacity shows a peak at some optimal 
power. This peak is at slightly higher power than that for 
the coherent capacity, because phase/timing fluctuations 
are more sensitive to XPM than amplitude fluctuations. 

It seems plausible that it should be possible to express 
the capacity of a channel purely in terms of its multiple 
order Green's functions. Unfortunately, we do not have 
such an expression. Indeed, only for a channel that is 
entirely determined by its two-point (or single particle) 
Green's function is it possible to write down such an ex- 
pression. In this case the channel is Gaussian and the re- 
sult is Eq.(28) as obtained by Shannon!!!. Obtaining such 
expressions would be a significant contribution to the un- 
derstanding of communication in non-linear media. The 
bound on coherent communication and the estimate of 
capacity for incoherent communication given above give 
an indication of the importance of including higher or- 
der Green's functions. Vertex corrections to higher or- 
der Green's functions determine the extent to which the 
channel has non-Gaussian statistics. Including these cor- 
rections can dramatically modify the functional depen- 
dence of ones estimates of system performance upon sys- 
tem parameters. 



VII. CONCLUSIONS 

In conclusion, we have considered the impediments to 
communication caused by non-linear interaction between 
the sub-bands of a WDM optical fiber system. These 
interactions between sub-bands lead to non-Gaussian in- 
put/output statistics. We have provided a simple char- 
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acterization of these non-Gaussian statistics in terms of 
vertex corrections to higher order Green's functions of 
the channel. Propagation of light in the WDM system 
under the action of XPM is described by a Schrodinger 
equation with a spatially and temporally random poten- 
tial. This potential encodes the effects of XPM. We have 
calculated the Green's function for propagation in this 
channel using Feynman path-integral and diagrammatic 
techniques. In the case where the signal correlations are 
approximated by delta-functions in time, these Green's 
functions may be obtained in closed form. The result 
shows a super-diffusive spreading of a signal pulse and 
of its arrival time. We have interpreted these results 
in terms of the amplitude and timing jitter of received 
pulses; both of which are standard measures of optical 
system performance. We find that the received signal 
pulse shows a super-diffusive spreading of its distribu- 
tion of arrival times; (t 2 ) cx L 3 . Finally, we have used 
the results to obtain lower bounds and estimates of the 
capacity of the WDM system. These are based upon the 
Shannon result for the capacity of a Gaussian channel. 
In the case of coherent communication, they lead to a 
strict lower bound upon the capacity. In the case of inco- 
herent communication, where only the signal amplitude 
and not its phase is used for communication, we obtain 
an estimate of the system capacity. The different de- 
pendences of our estimates of capacity for coherent and 
incoherent communication at large distances, exponen- 
tial and power-law, respectively, show the importance of 
accounting for the non-Gaussian nature of the channel. 



APPENDIX A: CALCULATION OF AMPLITUDE 
FLUCTUATIONS 

The results contained in Eqs.(22) may be derived using 
rules for the composition, 

J dt'g(xi,x';ti,t')g(x',x2- 1 t',t 2 ) = g(xi,x 2 ;h,t 2 ), 

(Al) 

unitarity, 



and time-reversal, 



G*(xi,x 2 ;ti,t 2 ) = Q{x 2 ,xi;t 2 ,ti) 



(A3) 



of the Green's function for a particular realization of the 
potential, V(x, t). These properties of the Green's func- 
tion may be confirmed by considering the path integral 
formulation, Eq.(6). 

Using these relations, the first of Eqs.(22) may be de- 
duced by the following series of manipulations: 



(E out (t)) = (A / dt'g(t,t') + n(t)) 



A( / dt'g(t,t')) 



= A(g(Lu = o)) 

This is nothing more than the average of Eq.(3). The 
second of Eqs.(22) requires a few more manipulations; 



(\Eout{t)\ 2 ) 



A 2 ( J dhdhgitMWitM)) + (Ht)\ 2 ) 
+A(n(t) J AiG*(t,ti))+A(n(t)* J dt 1 g{t,t 1 

a 2 J dtdhdt 2 g{t, h)g* {t, t 2 )/s{o) + N6(o) 
A 2 J dtdhdt 2 g(o, 0; t 2 , h)/S(o) + ns(o) 

A 2 J dtdhdt 2 5(ti - h)/6(0) + NS(0) 
(A 2 + N)S(0) 



We have used Eq.(A3) in moving from the first to second 
line, followed by Eqs.( Al ) and ( A2) moving between lines 



g(x = 0;t u t 2 ) = 6(h-t 2 ) 



(A2) 



two and three and three and four, respectively. We have 
also used the fact that the averages over the noise and 
over the potential are temporally invariant. Eq.(22c) may 
be deduced as follows: 



(5\E out (t)\ 2 S\E out (t)\ 2 ) = (\E out (t)\ 2 \E out (t)\ 2 ) - (\E out (t)\ 2 ) 2 

= A 4 J dt l dt 2 dt 3 dt A {g*{L-,tM)Q%^t2)Q{L]tM)g{L-tM)) 

+ANA 2 J dhdt 2 (g * (L;t,t 1 )g(L;t,t 2 )) + [2N 2 - (A 2 + N) 2 ] 5(0) 2 
« A 4 f dt 1 dt 2 dt 3 dt 4 (g* (L; t, ti)g(L; t, t 3 ))(g*(L; t, t 2 )g(L; t, U)) + [N(2A 2 + N) - A 4 ] 6(0) 2 



A 4 



dtdhdt 2 (g*(L; t, h)g{L- 1, t 3 ))/S(0) 



+ [N(2A 2 + N) — A 4 ] 5{0f 
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rj N(2A 2 + N)S(0) 2 



The third line is simply a rearrangement of the various 
terms arising after substitution of Eq.(3). A complete 
calculation of the first term in this expression requires 
evaluation of the 4-particle Green's function. We do not 
carry out this calculation here. Instead we make the 
approximation discussed in the text and ignore vertex 



corrections. This approximation is embodied in the step 
between the third and fourth lines. In moving to the 
fifth line, we have used the fact that the average over the 
noise and random potential leads to tempor ally in vari ant 
expressions. We may then use Eqs.(Al,[A2) and ([A3] ) in 
the same way as before in order to derive our final result. 
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